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У статті запропоновано метод вибору коефіцієнта масштабування недіадного  вейвлет- 
перетворення на основі аналізу частотного спектру вихідного сигналу. Розглянуто застосування методу 
на прикладах модельного та реального сигналів. 
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Вступ 

Вейвлет-перетворення (ВП) є особливим типом лінійного перетворення сигналів і 
даних, що ними відображаються, про процеси та властивості природних і штучних 
середовищ та об'єктів. За допомогою ВІ можна здійснювати частотно-часовий аналіз, 
результати якого містять як розподіл енергії сигналу за частотними складовими, так й 
інформацію про певні локальні координати, на яких відбувається швидка зміна 
відповідних складових сигналу. Наразі аналіз шляхом застосування ВП є одним із 
базових способів дослідження часових рядів, зокрема, й тих, що містять вихідні дані 
для систем підтримки прийняття рішень. 

Постановка проблеми 

На практиці найчастіше застосовуються діадні дискретні ВП, тобто ВП з коефіцієнтом 
масштабування, що дорівнює 2, бо вони мають найбільш ефективну програмну реалізацію. 

Проте в деяких випадках вибір коефіцієнта масштабування рівний 2 не є 
найкращим рішенням. Наприклад, якщо пов'язані особливості досліджуваного явища 
містяться одразу в кількох сусідніх частотних інтервалах діадного ВІЇ, то для 
адекватного аналізу цих даних набагато краще буде виконати розбиття частотної 
області на такі інтервали, які б повністю охоплювали зазначені особливості |2|. Інше 
міркування на користь недіадних ВП полягає в тому, що двійкове масштабування не 
завжди є природним для певної предметної області чи розв'язуваної задачі. 
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Оскільки замість двійки можна використовувати як коефіцієнт масштабування 
будь-яке інше раціональне число, більше одиниці, то постає задача оптимального 
підбору коефіцієнта масштабування. 

Основні поняття 

Нехай М -- раціональне число, а Ф(х) та м(х) -- функції з простору 12( 97). Функція 
ФО) називається функцією масштабування, якщо вона задовольняє співвідношенню (1): 


Ф(х) «М - У п Ф(Мх-к). а) 
кед 
Число М називається коефіцієнтом масштабування ВП, а набір коефіцієнтів (Лк) - 
низькочастотним фільтром розкладу. 
Вейвлет у(х) задовольняє співвідношенню (2), аналогічному до (1): 


м(х)-УМ- У яр 'м(Мк-К). (2) 
кед 
Набір коефіцієнтів (2) називається високочастотним фільтром розкладу. 
Якщо позначити через ак; та ак; відповідно апроксимуючі та деталізуючі 
коефіцієнти к-го рівня розкладу, то початковий сигнал 5(7) можна представити як суму 
його апроксимації на К-ому рівні 1 К деталізуючих складових: 


- з К - - 
«(а УМ задом Кк) У УМ 4 (міт), 
) | 

Аналіз останніх досліджень і публікацій 

На даний момент запропоновано декілька методів для виконання недіадного 
дискретного ВП із різним математичним підгрунтям та характеристиками. У роботі |3| 
авторами був проведений аналіз методів та визначені їх особливості. 

Критерієм якості розкладення пропонується обрати ентропію. Інформаційна 
ентропія - це міра невизначеності чи непередбачуваності інформації |4). У роботах |3| 
та |6| саме ентропія використовується як критерій для отримання оптимальної 
декомпозиції сигналу. 

У |71 ентропія вейвлет-розкладення У/Е визначається наступним чином: 


М 


МЕ а-, р) п ру, 


уні 
де М - максимальний рівень розкладення; 
рі - відносна енергія вейвлет-розкладення, яка показує розподіл енергії по рівнях розкладення: 


Е; - величина енергії на )-ому рівні розкладення з кількістю деталізуючих 
коефіцієнтів М. Сукупність величин Е; утворює спектр вейвлет-енергії вихідного 
сигналу: 
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Ео - повна енергія вейвлет-спектру: 


М 
ЕЕ 2 В, і 


зі 


Мета дослідження 

Мета дослідження полягає в розробці методу вибору на основі частотного спектру 
вихідного сигналу такого коефіцієнта масштабування ВП, при якому буде досягнуто 
найбільш точне виділення наявних у сигналі періодичних особливостей. 

Опис запропонованого методу та експериментальних результатів 

Вибір коефіцієнта масштабування пропонується виконувати на основі аналізу 
частотного спектру вихідного сигналу. Періодичним особливостям відповідають 
частоти у спектрі із явно вираженими максимумами амплітуди. Отже, коефіцієнт 
масштабування потрібно (за можливості) обирати таким чином, щоб до кожного 
частотного інтервалу деталізуючих складових на кожному рівні розкладення потрапляв 
тільки один з таких максимумів. 

У даній роботі недіадне ВП здійснюється методом, запропонованим в (8). Його 
основними перевагами є можливість використання як коефіцієнта масштабування 
довільного раціонального числа р/цд так і простого пірамідального алгоритму вейвлет- 
розкладення у частотній області. На кожному рівні розкладення обчислюються 
коефіцієнти для однієї апроксимаційної та р-д деталізуючих складових. За рахунок 
кількості деталізуючих складових та довжини відповідних частотних інтервалів можна 
досягти більш точного відокремлення особливостей вихідного сигналу. 

Найбільш перспективним для проведення подальшого аналізу будемо вважати 
таке значення коефіцієнта масштабування, при якому ентропія вейвлет-розкладення 
буде мінімальною. 

Розглянемо спочатку модельний сигнал (рис. 1), який представляє собою суму 
двох синусоїдальних сигналів. Його частотний спектр наведено на рис. 2, із якого 
видно, що обидві частотні особливості потрапляють в один інтервал діадного ВП. 

Тому при використанні для аналізу цього сигналу діадного ВП чітке 
відокремлення зазначених періодичних особливостей неможливе. 

Результати застосування вейвлетів Добеші 4-го порядку та Хаара наведені на рис. 
3 та 4 відповідно. Як можна бачити, жодна з отриманих при цьому деталізуючих 
складових не схожа на синусоїду. 
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Рис. Модельний сигнал 
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Рис. 2 Частотний спектр модельного сигналу та його діадне розбиття 
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Рис. 3 Застосування вейвлета Добеші до модельного сигналу 5: 


а -- апроксимаційна складова 5-го рівня розкладення, а), ..., д5 - деталізуючі 
складові відповідних рівнів розкладення 
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Рис. 4 Застосування вейвлета Хаара до модельного сигналу 5: 
а -- апроксимаційна складова 5-го рівня розкладення, а), ..., 45 - деталізуючі 
складові відповідних рівнів розкладення 


Проте, якщо виконувати розбиття частотного спектру на три частини (тобто 
виділяти при ВП на кожному рівні дві деталізуючі та одну апроксимаційну складові), 
то картина змінюється (рис. 5). Один з частотних піків потрапляє в частотний 
проміжок, що відповідає одній з деталізуючих складових першого рівня, а другий - в 
частотний проміжок для іншої з деталізуючих складових другого рівня вейвлет-розкладу. 

Виходячи з цього, застосуємо до модельного сигналу обране недіадне ВП із 
коефіцієнтом масштабування, рівним трьом (тобто покладемо р-3, 4-1). Обмежимося 
також третім рівнем розкладу (рис. 6). Як можна бачити, сума першої деталізуючої 
складової першого рівня та другої деталізуючої складової другого рівня розкладу 
достатньо близька за своєю формою до початкового модельного сигналу (рис. 7), 
відповідне середньоквадратичне відхилення дорівнює 0.1. Таке значення відхилення 
викликано впливом крайових ефектів перетворення Фур'є. Якщо відкинути по три 
відліки на початку та в кінці сигналу, то середньоквадратичне відхилення зменшується 
в 3,37 рази. 

Той факт, що недіадне ВП виявилося кращим за діадне, підтверджується і 
чисельним критерієм. Величина ентропії для тріадного розкладу дорівнює 0,6817, а для 
діадного - 0,9335 та 0,9051 для вейвлетів Хаара та Добеші четвертого порядку відповідно. 

Розглянемо тепер приклад з реальними даними, а саме: середньорічну ціну на 
пшеницю за період з 1500 по 1869 роки |9| (рис. 8). Відповідний частотний спектр 
наведено на рис. 9. 
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Виходячи з візуального аналізу частотного спектру, було обрано декілька значень 
коефіцієнта масштабування, при яких виконано вейвлет-розкладення вихідного сигналу 
до третього рівня та обчислено значення ентропії вейвлет-розкладення. Результати 
обчислень зведено до таблиці І. 
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Рис. 5 Тріадне розбиття частотного спектру модельного сигналу 
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Рис. 6 Недіадне ВП із параметрами р-3, 4-1, перший рядок - вихідний сигнал, 
наступні рядки - складові першого, другого та третього рівня розкладення (зліва 
направо - дві деталізуючі та апроксимаційна) 
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Рис. 7 Модельний сигнал (зверху) та сума двох деталізуючих складових (знизу) 
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Рис.9 Частотний спектр реального сигналу 
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Таблиця 1 - Значення ентропії для обраних коефіцієнтів масштабування 


Коефіцієнт масштабування Ентропія 
2 (Добеші 4-го порядку) 0,7270 
2 (Хаар) 0,9598 
9/4 (-2,25) 0,6950 
5/2 (-2,5) 0,6352 
ЗП 0,6624 


Найменше значення ентропії вейвлет-розкладення було досягнуто при коефіцієнті 
масштабування 5/2. І дійсно, при цьому значенні коефіцієнта певні деталізуючі 


складові, а саме: друга на другому та третьому рівнях розкладу вказують на наявність 
періодичних особливостей по всій довжині сигналу (рис. 10). 


Виходячи з частотних інтервалів обраних деталізуючих складових, можна 
підрахувати, що цим особливостям відповідають періоди в 24-30 років та 60-70 років. 
Відповідь на питання про те, що саме призвело до таких періодичностей, потребує 
окремого дослідження і буде залежати, скоріш за все, від різних супутніх факторів 
(кліматичні умови, історичні події тощо). 
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Рис. 10 Обрані деталізуючі складові розкладення при р-5, 4-2 

Висновки 

У роботі запропоновано здійснювати вибір коефіцієнту масштабування дискет- 
ного ВП для виділення періодичних особливостей у вихідному сигналі на основі його 
частотного спектру. 

Запропонований підхід було проілюстровано на модельному та реальному 
сигналі. Отримані результати показали, що використання недіадного ВП із визначеним 
коефіцієнтом виділяє наявні в сигналі періодичності, в тому числі й не очевидні. 

Представляється перспективною розробка методу, в якому вибір коефіцієнту 
масштабування буде виконуватися автоматично на основі аналізу спектру вихідного 
сигналу. При цьому потрібно також приділити увагу співвідношенню між кількістю 
деталізуючих складових та шириною частотних інтервалів, що їм відповідають. 
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КЕ5ОМЕ 

О. Срегіоу, У.У. Маїспукоу 

Мешоа Гог деїегтіпіпе дйабіоп Гасіог ої (бе поп-дуадіс уауеіеї (гап5їогт Базей 
оп Фгедипепсу зресігит 

Іп із агісіе, ацібог5 ргорозед а піефоай Баз5ед оп Ше апаіузіз ої Ше Їгедиепсу 
зресігит ої Ше іприк 5іспа! їог дегегтіпіпя Ше дайоп Касіог ої ре поп-Фуадіс ухауєієі 
сап5їогт угрісі уії аспіеує Ше тозі ассигаїе деїесйоп ої Ше регіодісаї Геашшгез Ша аге 
рге5епі їп Ше 5ієпа!. 

Гацоп Тасіог баз (о Бе 5еіесіва їп 5исП а мау аг еасі йедиепсу Бапа аї еуегу Ісусі 
ої ууамеїеі дДесотрозбійоп 5пошіа сопіаїп по поге їбап опе плахітит ої ігедиепсу апарішає. 

Бог поп-дуадіс ууаусеіеі: десотрозійоп, Ше тео іфаї иц5е Ше агріїгагу гайопаї 
питрбег а5 а Фашоп Гасіог Па5 Бееп 5еЇесіед. Епігору ої ууаусіеі десотрозбійоп у/аз спо8еп 
аз а сгікегіоп ої десотрозійоп дцашу. 

Арріїсайоп ої ргоро5ед теїрой у/а5 детоп5ігаїсд оп (о ехагаріе5 уп плодеі апа 
геа! Чака. ТБе ауегаєе аппиаі ргісе ої м/Пеаї омег Бе регіо тот 1500 (о 1869 уеаг5 ууа5 
аКеп аз а геа! 5ієпа!. Тре ргоро5ед теїроа аПоууеа ідепійутпе їп (Фі8 зієпа! бо регіодісаїз8 
сусівз їп 24-30 уеаг5 апа 60-70 уеагя8, муБісП ухеге пог уї51б1е їп ре огієїпа! зієпаї. 

Ії зеет5 рготізіпе, (0 іпаргоуе бе пецодй їп мубісп (пе Фйайоп Гасіог 5еРесійоп млі! Бе 
допе аціотайсаПу Бу апаїугіпе, Ше Їгедиепсу 5ресігит ої Ше огієїпа! 51єпа!. 
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